function [roc,ac]=rocSong(a,c)

%Pour chaque valeur de confiance, proba d'avoir cette confiance sachant
%qu'on est correct, en fonction de proba d'avoir cette confiance sachant
%qu'on est incorrect :

% p (con?dence = i | correct) =f( p (con?dence = i | incorrect) )


roc=[];
    for j=1:6
        nc(j)=sum(c==j);
        ac(j)=mean(a(c==j));
        roc(j,1:2) = [ sum(c==j & a==1)/sum(a==1) sum(c==j & a==0)/sum(a==0) ]; 
    end
    roc = [0 0; cumsum(roc)];
    X=roc(:,1);
    Y=roc(:,2);